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METHOD OF RELATING MODAL DAMPING TO LOCAL DAMPERS IN 
LUMPED-PARAMETER SYSTEMS 
by Harry J. Koenig and Daniel I. Drain 
Lewis Research Center 

SUMMARY 

The longitudinal dynamic characteristics of a linear structure may be approximated 
by a lumped-parameter model, which usually includes only the stiffness and mass distri- 
bution of the system. In general, the effect of damping is not included as a lumped 
parameter but is treated as an overall or modal characteristic superimposed on the re- 
sults of the model. This report investigates the relation between the modal damping 
characteristics and the individual section dampers of a lumped-parameter system that 
could model a launch vehicle for longitudinal motion. The analysis presented shows that, 
for the typical damping range of launch vehicles, an approximately linear relation exists 
between local and modal damping. This linearity allows the determination of either local 
or modal damping from known values of the other as well as a method of correcting test 
stand data for the support system damping effects. Also, it can furnish insight into the 
relative effectiveness of each local damper on the modal damping at each resonant fre- 
quency. Two examples of the use of the analysis are included herein as well as a 
FORTRAN IV computer program suitable for calculations of the damping relations in 
lumped models up to 10 sections in size. 


INTRODUCTION 

A commonly used method for the dynamic analysis of structures is the development 
of a linear spring-mass model of the structure. The values of the masses and spring 
constants are derived directly from the weight and stiffness of the lumped section of the 
structure. The approximate resonant frequencies of the structure can then be found by 
solving for the resonant frequencies of the lumped parameter model. For a complete 
dynamic description of a structure, the damping characteristics must be known as well 
as the resonant frequencies. In the design stages or when the structure cannot be shaken 



experimentally, a modal damping is assumed through previous experience with similar 
structures. 

Although the modal damping is a valuable entity, it is not sufficient to predict the 
effect of a modification to the structure. Through the use of individual section damping, 
the change in the structure dynamics because of a modification would be predicted more 
easily. A second valuable use would be in modeling the structure on the analog computer 
for longitudinal stability analysis. 

This study investigated the relation between the individual section dampers and the 
modal damping of a series mass-spring system. The system considered would model a 
space vehicle structure for longitudinal motion. A method of solving for the damping 
characteristics of discrete systems has been proposed by DaDeppo (ref. 1) for the special 
case of damping proportional to either the mass or the spring matrix of the system. This 
report differs in that it is not limited to cases where special relations exist between the 
system constants, and it is written in terms of the physical coordinates of the system. 

The relation of each individual section damper to the modal damping was found in 
terms of damping influence coefficients at each resonance for lightly damped (oscillatory) 
systems. Two examples of this method are given in the section APPLICATIONS. One 
example illustrates how to compensate for the influence of the suspension system damping 
on the experimental test data from a shaker facility. 


ANALYSIS 

The development of the relation between modal and section damping will be presented 
in steps. The equation of motion for a generalized single mass will be developed first; 
then these equations will be applied to the multimass system to form the system equations. 
Following this will be the procedure to solve the system equations, a brief discussion of 
modal modeling, and then the matching of the system equations to the 
modal equivalent model. This procedure yields equations that contain 
the information necessary to establish the general relation between 
modal and section damping. The general equations will then be applied 
to the physical case of a grounded two-mass system. 

General Equations of Motion 

The coordinate system and the related elements of mass, springs, 
and dampers for the i^ 1 mass of a multimass system are illustrated in 
figure 1. The equation of motion about the i th coordinate may be 
written immediately as 


r* X: 



*- x i+l 


Figure 1. - Generalized 
coordinate system. 
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Q. = M.x. + B.(i. - x i+1 ) + K.(x. - x. +1 ) 




Collecting common coordinates and taking 
the Laplace transformation with zero initial 
conditions yield 

[m 4 S 2 + (B i _ 1 + B.)S + K._ 1 + Kjjx.(S) 

- (B.s + iq)x i+1 (s) 

Figure 2. - Boundary conditions. 

- (B. jS + K._ 1 )x._ 1 (S) = Q.(S) (2) 

th 

which is the desired form of the general equation about the i mass. 

In addition to the general equation, two boundary conditions consistent with a series - 
type system were considered. The two conditions are the special constraint about the 
top and bottom mass of the system, as indicated in figure 2. The equation about the top 
mass may be derived from equation (2) by using the constraints B^_ ^ = Kj_^ = 0; that is, 

(MjS 2 + BjS + K 1 )x 1 (S) - (BjS + K 1 )x 2 (S) = Q X (S) (3) 

The equation about the bottom mass may be obtained by using the constraint x ]i+ -^ = 0: 



(a) Top mass configuration. (b) Bottom mass configuration. 


(M n S 2 ♦ (B n ♦ B n . 1 )S + K,, + K n _j)x n (S) - (B n _,S + K^^fS) = Q„(S) (4) 

Application of General Equation of Motion to n-Mass System 

The equation of motion for an n-mass system, as illustrated in figure 3, was written 
in matrix form. For the problem considered, all generalized forces Ql were set equal 
to zero except Q n , the force on the bottom mass, where F(S) was substituted for Q^S). 
The application of a sinusoidal force to the bottom mass is analogous to longitudinal 
shaking of a launch vehicle in a test facility. The equation of motion of the n-mass 
system in matrix form is 
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Figure 3. - Typical 
n-mass system. 



From this matrix, the relation or transfer function between the 
force at the bottom mass and the motion of the top mass Xj(S)/F(S) 
was selected for investigation because of the absence of antiresonances 
or node points. From equation (5), Xj(S)/F(S) may be solved for by 
Cramer's rule (ref. 2), and the resulting solution will have the general 
form 


4 







( 6 ) 


x^S) 

~¥&) 


(B^ + K 1 )(B 2 S + K 2 ) . . . (B n _ 1 S + y 
a 2n s2n + a 2n-i s2n_1 + ■ * • + a l S + a 0 


In the general form of the transfer function, the coefficients in the numerator are 
the section dampers and spring constants. In the denominator, or characteristic poly- 
nomial, the coefficients a 2n , a 2 n _i> • • • > a o are functions of the springs, masses, and 
dampers. A method of evaluating these coefficients is provided in appendix D. 


Modal Modeling 

The classical approach to modal modeling is to solve for a set of independent modal 
solutions, the sum of which is equal to the transfer function of the system. In order to 
do this, a coordinate transformation must be made. This transformation is necessary to 
obtain the independence or orthogonality of the solutions; however, it limits the useful- 
ness of the classical approach. 

In order to avoid these limitations, the system transfer function was matched at each 
resonate peak by a second-order spring -mass -damper system, as illustrated in figure 4 


Lumped-parameter 
model of actual 
system 



Figure 4. - Relation of lumped-parameter model to its modal equivalent models. 
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for a five-mass system. The mass, damper, and spring of each equivalent-system are 
represented by M , B , and K respectively; the additional Roman numeral sub- 

•L «L 

script indicates the order of the mode. Thus, the modal transfer function evaluated at 
the first resonant mode o>j would be 


Cl)t M tS + B jS + K j 

I eq, r eq, Jr eq, I 


x(S) 

F(S) 


To complete the discussion of modal modeling, it is necessary to relate the modal 
constants with the frequency response characteristics of the system. This can be done 
readily when the damping ratio | of the system is on the order of 10 percent or less, 
and a plot of the amplitude as a function of frequency characteristics of the system has 
been obtained. Since the modeling is done at resonant condition, the resonant amplitude 
ratio and frequency a> n are related to the modal equivalent values of mass, spring, and 
damper by the following equations: 


B 

M 


eq _ 


= 2 &n 


eq 


i oj=0 


u 


D M eq 


lu)=0 


K 


eq 


Combining equations (8) to (11) yields 





( 8 ) 

(9) 

( 10 ) 


( 11 ) 


( 12 ) 


These are approximate equations for a second-order system; the assumption used to 
obtain these equations is £ < 0. 1 and, therefore, the peak, damped, and undamped 
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frequencies are very nearly equal. 

With the aforementioned definitions and equations, the relation between the section 
dampers and modal damping can be developed. The relation will be developed analytically 
first and then applied in two examples. 


General Damping Equation 


In the development of the final equation of this section (eq. (18)), three approxima- 
tions will be used, the justifications of which will be discussed in the section Linearity 
Assumptions. The equation development is started by rewriting equation (6): 


Xj(S) 

F(S) 


(BjS + fqXBjjS + K 2 ) . . . (B n _ 1 S + K n _ x ) 


*2n s2n + *2rt-l S 


2n-l ^ a 2n-2 , 

+ a 2n-2 S + • 


( 6 ) 


+ a^S + a Q 


Substituting jo> for S, collecting the real and imaginary terms, and taking the magnitude 
yield 


Xi(jw)| 


F(ju>) 


(jBjU) + K 1 )(jB 2 w + K 2 ) . . . (jB n _ 1 w + K n _ 1 ) 


2 n 2 2 

_ a 2n( -w ) +a 2n-2^ -w ^ +- ' ' + a 0_ "^L^n-l^ ) + - • • +a l 


c o 


(13) 


,th 


Evaluating at which is the peak frequency of the fi mode, permits the use of the 
approximation that the real portion of the denominator of equation (13) at a resonate peak 
is zero; that is, 


/ 2 \ n / 2 \ n- l 2 

a 2nV _a W + a 2n-2\ _&, S2/ +. • . + a 2 (-o> ) + a Q = 0 

Substitution of equation (14) into equation (13) results in 


(14) 


x iG“> n ) 


^ B 1 w S 2 + K l)(3 B 2 a, J2 + K 2 ) • • ' 

(3 B n-l w O + Vl> 

F <J<o 0 > 


j[ a 2n-l (" w n) + a 2n-3 (“4) + • 

. . +a 3 (-o4) + a x 


(15) 


2 / \ 2 

Assuming B? « and inverting equation (15) yield 
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( 16 ) 



However, the left side of equation (16) is, from equation (12), equal to B eq 
tion (16) is, therefore, the general damping equation; that is, 


1 

/ 2 \ n_1 

( 2 \ n “ 2 

B eq,fi n _i 

a 2n-l\ _a W + a 2n-3 ' 

+ • • * + a l 


i=l K ‘ 


Equa- 

(17) 


In equation (17), the coefficients a. 2 n ^, . . . , a^ are linear functions of the local 
dampers B. if it is assumed that in each a term the sum of the terms in Bj is much 
greater than the sum of the product terms (e.g. , B^Bj ) . With this assumption and a com- 
bination of like terms, equation (17) becomes a linear form that yields 

B eq,n = K B 1 + a 2 B 2 + • • • + <*i B i + • • • + «n B nl < 18 > 

where a n is a combination of the systems masses, springs, and resonant frequencies. 


Linearity Assumptions 

In order to obtain the linear damping equations, three approximations are made. All 
approximations are consistent with light damping, as will be shown in this section. The 
first approximation, at the peak frequency a>^, is 


l 2n 


(-<4) + a 2n-2 ( -w n) 



(14) 


The validity of this assumption can be argued from the limiting case of no damping. When 
there is no damping, the resonances are infinite and there will be no imaginary portion 
of the denominator in equation (13). Thus, the denominator consists only of equation (14), 
and it must go to zero in order to get an infinite resonance. The addition of a small 
amount of damping has two effects: the imaginary portion of the denominator in equa- 
tion (13) exists and limits the resonant peak of the equation, and there is a slight modifi- 
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cation in the a coefficients of equation (14) (appendix B), which will cause a small shift 
in the resonant frequency at which equation (14) goes to zero. This shift in resonant 
frequency is assumed to be negligibly small and is valid only for small values of damping. 
The second approximation is 


4 


— i = 1, 2, 3, . . . , n - 1 


(19) 


The second approximation is used to eliminate the section damping from the numerator of 
equation (6). The effect of this assumption is to replace the product terms of B^S + 
by a product of K. 's in the calculations of B e ^ and M e q. This relation generally holds 
for the lower resonant frequencies; it may not hold for the higher resonant frequencies, 
however, thus introducing errors into the solutions. 

The third approximation is that each a coefficient in the imaginary polynomial 
(eq. (17)) is a linear function of the section dampers. This approximation is valid, as 
will be shown by example, if 

^ Terms linear in B^ » Product of B^ terms (20) 


where 


i = 1, 2, 3, . . . , n 

This approximation was investigated for systems of up to seven masses. The terms that 
were linear in B outnumbered greatly the product of B terms. For light damping 
(i.e., small values of B), the products of the B terms are second- and third-order 
effects, which add to the validity of the approximation. As an example in a three-mass 

O 

system, a product of B terms appears as a coefficient of S (see appendix C). The total 
coefficient is 

a 3 = (MjK 2 + M 1 K 3 + M 3 K 2 + M 2 K 2 + M 2 K 3 )B 1 + (M 1 K 1 + MjKg + Mgiq + M 2 K 3 )B 2 

+ (M 1 K 1 + + M 2 K 1 )B 3 + BjBgBg (21) 

In this coefficient there are 12 terms that are linear in B (assuming that M and K are 
constant). There is one product of B terms, BjBgBg. For light damping, the product 
of any mass spring MK is much larger than the product of damping coefficients BB. 


9 



The two other coefficients to be used in equation (17) for a three-mass system are 


a^ = MjMgBg + M^MgBg + M-jMgB^ + MjMgBg + 


( 22 ) 


a l = k 2 k 3 b i + Kl K 2 B 3 + Kl K 3 B 2 


(23) 


Coefficients a^ and a^ are linear in the section dampers without approximation (20). 
With the heuristic argument as given for coefficient a 3 , the validity of the third approxi- 
mation may be carried to an n-mass system. 

Each of the approximations has been discussed in general to indicate how each will 
affect the linearity of equation (18). A specific example, however, shows the combined 
effects. In this example, the section dampers were all made equal, and the normalized 

B for each mode was calculated for increasing values of section damping. The system 
eq 

model was the five-mass system used in the APPLICATIONS section of this report. The 
results are illustrated in figure 5 and indicate four distinct characteristics: 

(1) The straight-line slope at low damping is identical to that predicted by the lin- 
earized equation (18). 

(2) The first mode slope (rigid body) is a constant. 

(3) The second and fourth mode slopes are lower than the initial slope for large 
values of damping. 
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(4) The third and fifth mode slopes rise above the initial slope for large damping. 

The variances in slopes can be explained as follows: the rigid body mode behaves 
like a standard second-order system and, thus, should be linear; the sagging slope shows 
the influence of the violation of the second approximation; and the rising slope results 
from the resonant frequency shift due to damping . In this example, the frequency shift 
of the second and fourth modes is about 1 percent, while that for the third and fifth modes 
is 10 and 20 percent, respectively, at a damping value of 4000 pounds-second per foot. 
Thus, the two effects can be separated. 

In figure 5, the normalized B can be considered the critical damping factor of 
a second-order system. Thus, for this example, the linearized equation concept is rea- 
sonably accurate for all modes up to 5 percent modal damping and, except for the fifth 
mode, good to 10 percent of critical damping. 


Damping Influence Coefficients 


Around the resonances of a multiordered system, the modal damping can be approxi- 
mated by a linear function of the section dampers (Linearity Assumptions). With the use 
of this linear function, the idea of damping influence coefficients may be introduced. 

The damping equation about the resonance S2 is 


«1 B 1 + “2 B 2 + 


“i B i + 


+ 0 
n n eq. 


(18) 


where (i = 1, 2, . . . . , n) are constants and called damping influence coefficients. 
Due to physical consideration, the of equation (18) must be of the same sign or 
sustained oscillations could exist in a dissipative system. Therefore, may be 
defined more explicitly as 


1 B i 


(24) 


where 


— 0 k — 1, 2, 3, . . . , n*. k ^ i 

By using this process (i. e. , setting all values of B equal to zero except the B of in- 
terest), the damping influence coefficients may be obtained with a digital computer. The 
programs used are given in appendix D. The influence coefficients are combinations of 
the systems masses, springs, and modal frequencies, as will be shown in the following 
sections. 
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Development of Damping Equations for Two-Mass System 


For a two-mass system, the coefficients of the characteristic polynomial are found 
by taking the determinant of the system matrix of equation (5). The matrix for the two- 
mass system is 


MjS> + BjS + Kj 


-(BjS + Kj) 


-(BjS + K x ) MgD + (Bj + B 2 )S + K x + Kg 


(25) 


In order to obtain the characteristic polynomial, the determinant of equation (25) is 
found. The polynomial generated is of the form of the denominator of equation (6), that 
is, 


a 4 S 4 + a 3 S 3 + a 2 S 2 + a x S 


where 


a 4 = M 1 M 2 

(26) 

a 3 — - M j + M j B + M ^ )1 - 

(27) 

= Mj^Kj + MjKg + M 2 K x + BjBg 

(28) 

a l = K 2 B 1 + K 1 B 2 

(29) 

a 0 = K 1 K 2 

(30) 


The linear damping equations for a two -mass system may be obtained from the 
general equation (eq. (17)): 


B an t = — 

eq , I Ki 

B - 1 
B eq,n-- 



(31) 


( 32 ) 
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When equations (27) and (29) are combined with equations (31) and (32), the following are 
obtained: 


B M T = — 
eq, I 


eq, n 


K, 


[k 2 - (Mj + ♦ (K, - Mj^Db. 

[k 2 - (Mj + M 2 )4jBj + (Kj - M 1 <4 )b. 


(33) 


(34) 


From equation (33), it can be seen that the damping influence coefficients for the first 
mode are 


a l = 


K, 


[k 2 - (Mj + M 2 )o|] 


(35) 


a 2 


K, 


( K i - M i“i 2 ) 


(36) 


The physical significance of the influence coefficient is that is the amount of modal 
equivalent damping contributed by Bj when Bj has a unit value of damping (i.e., 

B = 1 lb-sec/ft). 


APPLICATIONS 

After the linear damping equations and modal dampings are found, a system of n 
simultaneous nonhomogeneous linear equations is obtained. These equations have the sec- 
tion dampers as their variables and may be solved simultaneously. 

The solution to the system of equations is the actual value of section dampers, which, 
if substituted into the spring -mass model, would give a reproduction of the amplitude 
frequency plot. A Gaussian reduction program is supplied in appendix D for the solution 
of the simultaneous equations. Two examples of the general technique are given in the 
following sections. 


Example 1 

K the spring -mass system (fig. 6) with the parameter values of 
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= 1 


M 2 = 2 


Hi 



K x = 20 K 2 = 10 


B 1 = 0. 1 B 2 = 0. 05 


is given, the plot of amplitude as a function of frequency may 
be calculated directly from the theoretical transfer function of the 
system. This amplitude -frequency plot is presented in figure 7. 

For this example, it could be assumed that the plot is data and Bj 
Figure 6. - Two-mass system, and B 2 of the aforementioned model are unknown. The problem 

is finding the section dampers, given the spring -mass model and the 
amplitude frequency plot. 

The two linear damping equations for the fourth-order system are 




Kh 


[k 2 - (Mj + M 2 )o| B x + (Kj - M iW j)b, 


B eq,n - [ k 2 - <“i + M 2>“n] B ! - ( K 1 - M i"n) B ; 


(33) 


(34) 


Four values must now be obtained from the experimental amplitude -frequency plot 
(fig. 7). These values are the peak amplitudes and resonant frequencies: 

|x |p^ j = 12. 5 coj = 1. 77 



[x |p n = 0. 385 Wjj = 5. 64 


Substituting the values of 
results in 

Kp K 2 , Mp M 2 , o>j, and Wjj into equations (33) and (34) 



B._ T = |0. 0292B, + 0. 8430B 9 1 
eq, II l & \ 

(37) 


B eq,n- I-4.279BJ- 0.6980B2I 

(38) 
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20 1 


10! 

8 

6 



nl | | | \ I | I . I 1 1 I ' — » « > 

o 1 2 3 4 5 6 7 

Angular frequency, rad/sec 

Figure 7, - Amplitude-frequency plot for fourth -order system. Parameters: 
mass, - 1 and M 2 * 2; spring constant, Kj - 20 and K 2 - 10; damping, 

* 0. 1 and B 2 * 0.05. 

From equation (12), 


eq, 1 (1. 77)(12. 5) 


O. 0451 


(39) 


eq ’ n (5. 64)(0. 385) 


= O. 461 


(40) 


Solving equations (37) and (38) yields 

Bj = 0.102 B 2 = 0.050 


which are the values of section dampers that would generate figure 7 with the given 
spring -mass model. 

As the number of sections increases, the problem becomes more difficult, and so- 
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lution by hand, as in example 1, must be abandoned. Example 2 has five sections, and 
the problem is solved by using the digital computer. 


Example 2 


Part A . - A problem that might be presented in practice would be the evaluation of 
the section dampers of a launch vehicle if the following conditions had been filled: 

(1) The vehicle had been shaken in a test stand. 

(2) The plot of longitudinal amplitude as a function of frequency had been obtained 
(fig. 8). 

(3) The spring -mass model had been obtained from the weight and spring distribution 
of the vehicle. 


TABLE I. - ASSUMED FIVE-MASS 
SYSTEM PARAMETERS 


Section 

Mass, 

slug 

Spring, 

lb/ft 

Damper 

lb-sec/fl 

1 

1060 

8. 65X10 6 

500 

2 

68 

5.35 

900 

3 

58 

14.0 

640 

4 

98 

18.0 

1300 

5 

178 

.21 

1600 


In order to demonstrate the technique to be used 
in the solution to this problem, figure 8 was calculated 
from the transfer function of a five -mass system with 
the parameters given in table I (see fig. 4, p. 5) and 
was assumed to be the experimental data. Presuming 
that the experimental data in figure 8 and the spring- 
mass model are available, the problem is to find the 
dampers in the vehicle that are characteristic to the 
generated plot. 
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The main program (appendix D) is used, and the lumped-mass and spring constants 
are applied to develop the following linear damping equations: 

B e q j = 0. 0003 IB x + 0. 00089B 2 + 0. OOOI6B3 + 0. 00011B 4 + 0. 93354B 5 ' 

B eq H = °* 52103B l + !• 2987B 2 + 0. 13964B 3 + 0. 03485B 4 + 3. 5894B 5 
B eq m = 47 • 788B 1 + 2 - 1330B 2 + 14. 120B 3 + 16. 146B 4 + 7. 2588B & ^ (41) 

B eq IV = 93 * 375B 1 + 256 * 54B 2 + 1 * 4449b 3 + 77 • 778 B 4 + 11. 395B & 

B eq y = 87. 102BJ + 2389. 2B 2 + 5087. 6B 3 + 1507. 7B 4 + 51. 669B & 

Inversion of the matrices of equation (41) results in 

B 1 = -2. 732B eq> j + 0. 6818B eq> n + 0. 01937B eq> m - 0. 00293B eq> ^ - 0. 715xl0"*B eq> y 
B 2 = -1. 923B eq> j + 0. 5123B eq> n - 0. 008136^ m + 0. 001302B eq> jy + 0. 819xl0“ 5 B eq> y 
B 3 = -1. 876B eq> j + 0. 4938B eq> n + 0. 002445B eq> ffl - 0. 00416B eq> Iy + 1. 77 xlO“ 4 B eq y 
B 4 = 9. 502B eq> j - 2. 5176B eq> n + 0. 003516B eq> m + 0. 01216B eq> jy + 0. 556xl<T 4 B eq> y 

B 5 = 1 * 073B eq, I * °* 000501B eq, n + 4x10 ^eq, m “ 01x10 5B eq, IV " °* 191x10 7 B e q, V 


\-4 t 


(42) 


The experimental values of B__ must be found now from the experimental amplitude - 
frequency plot (fig. 8). From this plot, the resonant amplitudes and frequencies were 
obtained and are listed in table II along with the values of B eq ^ calculated from equa- 
tion (12). 

The equivalent dampings are then substituted into the inverted matrix, and the sec- 
tion dampers are calculated. In table m, the calculated damping values are compared 
with the assumed damping values of table I. 

A perfect agreement between the calculated and assumed section damper values 
should not be expected, since there were errors in reading the values of B gq from 
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TABLE n. - SYSTEM MODAL CHARACTERISTICS 


TABLE HI. - COMPARISON OF 


BASED ON FIGURE 8 


CALCULATED AND ASSUMED 


Quadratic 

Resonant 

frequency, 

I 

11.6 

n 

105 

m 

390 

IV 

512 

V 

734 


Amplitude 

Second-order 
modal system 
damping, 

B eq 

10 ”4’ 2 4 

1495 

IQ-5. 89 

7390 

10-7.42 

67500 

10-8.29 

381000 

10-9.7! 

6990000 


DAMPING VALUES 


Section 

damper 

Calculated 
damping 
using linear 
damping 
equations 

Assumed 
section 
damper 
used for 
figure 8 

B i 

641 

500 

B 2 

916 

900 

B 3 

661 

640 

B 4 

860 

1300 

B 5 

1600 

1600 


TABLE IV. - EVALUATION OF CALCULATION AND READING ERRORS 


Case 

Error in reading second-order 
equivalent mode damping, percent 

Resulting error in damping, percent 


B eq, I 

B eq,n 

B eq,m 

B eq, IV 

B eq,V 

B i 

B 2 

B 3 

B 4 

B 5 

1 

0 

0 

0 

0 

0 

8.4 

2.2 

4.4 

6.7 

0 

2 

. 1 

1.2 

-.3 

-2.2 

-2.4 

28.2 

1. 8 

3.2 

33. 8 

0 


figure 8 and since the analysis uses approximations. Two cases of the expected errors 
for this example are cited in table IV. Case 1 indicates the analysis approximation 
errors, and case 2 is the data used in this example. The large errors in the values of 
B^ and for the second case are not as important as might be expected. In fact, they 
indicate the lack of influence of these dampers on the transfer function. This lack of in- 
fluence can be demonstrated by calculating the transfer function with the results of this 
second example and comparing it with the original function. When this is done, the 
curve shapes are identical. The amplitude of the peaks for modes I to V are in error by 
0. 0, 1. 1, 0. 1, 4. 2, and 6. 7 percent, which is not much more in error than the original 
reading errors of the example. 

The masses and spring constants of example 2 are typical of an actual two-stage 
vehicle with conditions close to first-stage cutoff and supported in a shake stand. The 
symbols Bg and Kg represent the shake-stand suspension system damping and spring 
constant. Inspection of equation (41), which relates the modal to section damping, will 
show that the stand damping Bg constitutes 100 percent of the rigid body -mode damping 
Bgq j, as it should, but it also contributes 80 percent of the measured first-mode damp- 
ing B e q jj, which is undesirable. A reduction in the magnitude of the stand damping 
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from Bg = 1600 to Bg = 160 will reduce the stand damping contribution to the first mode 
to 25 percent, which is still substantial and which readily points out the need for com- 
pensating for the suspension system damping as the damping characteristics are being 
evaluated. Node point suspension, of course, is the ideal way to eliminate the effect of 
stand damping, but this, in general, is not practical for longitudinal dynamic tests. A 
method to eliminate the stand damping based on the technique presented herein is given 
in the next section. 

Part B . - By simply subtracting the stand damping, the spring -mass model may be 
used to extrapolate to the value of damping that would exist in a vehicle in free flight. 

The first step is to perform an analysis, as in Part A, which uses the B found from 

“ 4 . 

the amplitude -frequency plot. Equation (43) may then be written as 


B eq H ‘ 3 ' 5894B 5 = °- 52103B 1 + 1. 2987B 2 + 0. 13964B 3 + 0. 03485B 4 
B eq m “ 7 * 3588B 5 = 47 ■ 788B i + 2- 1330B 2 + 14. 12Bg + 16. 146B 4 


B eq, IV 
B eq, V 




11. 395Bg = 93. 375BJ + 256. 54B 2 + 1. 4449Bg + 77. 775B 4 
51. 669Bg = 87. 102BJ + 2389. 2B 2 + 5087. 6Bg + 1507. 7B 4 


(43) 


The linear damping equations, however, may be obtained for a free-free system as well 
as a grounded system. The method is the same as that used in appendix B, but the bottom 
spring and damper are set at zero. If the damping equations of the free-free system in 
Part A are generated, the linear damping equations obtained are 


B eq H ff = °- 4963B 1 + !• 248B 2 + 0. 13748Bg + 0. 03731B 4 
B eq m ff = 47 • 811B i + 2 - 069B 2 + 14. 124Bg + 16. 2808B 4 


> 


B e q iv ff = 93.908BJ + 257.31B 2 + 1.4069Bg + 77.819B 4 
B eq V ff = 87 • 349B 1 + 2395. 5B 2 + 5099. 4Bg + 1509. 3B4 

The similarity between equations (43) and (44) leads to the approximations 


(44) 
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( 45 ) 


B eq,n,ff - B eq,n " 
B eq,m,ff - B eq,m 


3.5894B 5 ^ 
- 7.3588B & 




B eq, IV, ff ~ B eq,IV ' 
B eq,V,ff = B eq,V ' 


11.395B 5 

51. 669B,- 

o J 


These equations illustrate that, even though the vehicle section dampers were not found 
with accuracy, the model of the vehicle will give the amplitude-frequency plot of the ve- 
hicle in free flight if the stand damping Bg is known. This is not an extremely stringent 
condition since Bg may be determined from the rigid body mode B e ^ j (eq. (42)). 

The influence coefficients can also be used as a guide to the most effective location 
for section damping of a vehicle during free flight. Thus, for the first free -free mode 
(eq. (44)), it can be seen that, by increasing Bg, the overall damping of that mode is in- 
creased optimally by 1. 248 times the damping introduced at Bg. If damping could be 
added to any section with equal facility (not true in practice), the optimum placement to 
increase overall effective damping would be at section 2. 


CONCLUSIONS 

For longitudinal oscillations of lightly damped structures, which are usually charac- 
teristic of launch vehicles, the relation of the overall modal damping to section damping 
has been shown to be approximately a linear function. A general procedure is furnished 
that will develop the influence coefficient of each section damper on the system damping. 
These influence coefficients thus can be used in an explanation of the effect of changes 
in local damping magnitude on the overall vehicle damping. They can also be used to 
provide insight into how the test vehicle suspension system damping will affect the meas- 
urement of vehicle damping characteristics in a test stand. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 21, 1966, 

180-31-01-01-22. 
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APPENDIX A 


SYMBOLS 


[A] 

matrix of order 2n 

« 

damping ratio, dimensionless 

a 

characteristic polynomial coeffi- 

P 

defined by eq. (D7) 


cient (see eq. (6)) 

<Pl 

first state variable defined as x, 

B 

damping, lb-sec/ft 

(n x 1) column vector 

[C] 

manipulative matrix (see 

<P 2 

second state variable defined as 

d 

appendix D) 

matrix component (eq. (5)) 

m 

dx/dt, (n x 1) column vector 
state vector, (2n x 1) column 

F 

force, lb 


vector 

P] 

identity matrix 

<j) 

angular frequency, rad/sec 

j 

yj-i 

Cl) 

n 

resonant angular frequency, 

K 

M 

n 

spring constant, lb/ft 
mass, slugs 

number of masses in system 

rad/ sec 
Subscripts: 

eq equivalent second-order (modal) 

system value 

Q 

generalized force, lb 

ff 

free-free system 

s 

Laplace operator, sec" 

i 

index of lumped-parameter 

X 

generalized coordinate, ft, 


system 1, 2, 3, . . . , n 


n x 1 column vector 

k 

dummy indexing variable 

a 

damping influence coefficient, 

J i 

index for y and exponent of [A] 

y 

M 

dimensionless 
defined by eq. (D8) 
number of resonances in system 

P 

O 

peak 

index of quadratic term (eq. (20)) 

i, n, m, . . ., fx 
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APPENDIX B 


DERIVATION OF MODAL EQUIVALENT MASS 


Grouping the odd and even powers of S from equation (6) yields 


XjCS) 

W) 


/a + a g2n-2 

\ a 2m + a 2n-2 & + 


(BjS + K 1 )(B 2 S + K 2 ) . . . (B n l S + K n l ) 
0 


+ <S 


2n-2 


2n- T" 


+ a 2n-3® 2n 4 + 


• + a l) S 

(Bl) 


The term 


a 2 n s2 " + a 2 n -2 s2 "' 2 + a 2 n -4 s2 ”' 4 + • • • + * a 0 < B2 > 

may be factored into the form 

a 2 n (® 2 + “f )( s2 + < 4 ) (® 2 + • • • ( S 2 + "m ) ® 3) 

where u>j, Wjj, are the frequencies where the real polynomial (eq. (14)) 

goes to zero. This may be shown readily when is substituted for S. With equation 
(B3), equation (Bl) may be rewritten as 


x 1 (S) 
F(S) ' 


a 2n(® 2 + 4 ) ( s2 + “m) ' ' ■ ( s2 + 4) S 2 + ( a 2n-l® 


<V + Kl)(B 2 S + K 2 ). ■ ■ (B n-1 S + K n . 1 ) 

■ + 3-^ S + 3 


2n-l 


2n 


(s 2 + 4) . . . (s 2 + «{)«s 


(B4) 


Equation (B4) is then conformed to equation (7): 

a 2 n ( s2 + " n )( s2 + " m )- • • ( s2 + a)2 ) 


M 


eq, I 


(B X S + K 1 )(B 2 S + K 2 )(B 3 S + K s ) . . . (B n _ 1 S + i^ j) 


S=i 


]Wj 


B 


eq, I 


a 2n-l s2n 2 + a 2n-3® 2n 4 + * * * + a l 


(B^ + K 1 )(B 2 S + K 2 ) . . . (B n l S + K^) 


S=j 




K eq,I = M eq,I W I 


(B5) 


(B6) 

(B7) 
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APPENDIX C 


COEFFICIENTS OF CHARACTERISTIC EQUATION OF THREE -MASS 

GROUNDED SYSTEM 

The coefficients of the characteristic equation of a three-mass grounded system are 
as follows: 


a 6 _ M 1 M 2 M 3 


(Cl) 


a 5 = MjMgBg + MjMgBg + MjMgB^ + MjMgB 2 + MgMgBj 


(C2) 


a 4 = MjMgKg + ^^MgKg + MjMgK^ + MjMgKg + + ^1®1®3 + ®^1®2®3 + ®1®2^2 

+ ^^2^3 ^ B2B3M2 + K^M2 Mj (C 3) 

3-0 = ^4^KgB 2 + M^K^B ^ + 3\4^K 2 Bg ^4gK 2 Bj. ^^2^2 B 1 

+ + B^BgBg + MgK^B 2 + MgK^Bg ^2^1^3 (C4) 

a 2 = M 1 K 1 K 2 + MjKjKg + M^KgKg + K^Bg + KgB 1 B 2 + MgKjKg + M^Kg + MgK^Kg 

+ K^Bg (C 5 ) 


a l = K 2 K 3 B 1 + K 1 K 2 B 3 + K 1 K 3 B 2 


(C6) 


a 0 ~ K 1 K 2 K 3 


(C7) 
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APPENDIX D 


DESCRIPTION AND LISTING OF DAMPER 

The main FORTRAN IV program and flow chart used in calculating the influence co- 
efficients are presented herein. The main program, DAMPER, has three subroutines 
called COEFF, GAUSSR, and POLYRT. Subroutines COEFF, which finds the coefficients 
of the characteristic polynomial, and GAUSSR, a Gaussian reduction subroutine, are 
presented herein. The subroutine POLYRT, which is a common polynomial root finding 
subroutine, is not supplied. 

The technique used in generating the influence coefficients is to set all dampers of 
the system to zero except the one in question, as is stated in the text in the sections 
General Damping Equation and Damping Influence Coefficients; the corresponding B 

eq 

is then calculated. 

The flow chart for DAMPER is given in figure 9. The FORTRAN IV program for 
DAMPER is presented after figure 9. 


Start 



Figure 9. - Main program-DAMPER. 
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o o n n n 


01 HENS ION QOC20) ,V( 11) ,6B( 10,11) ,B( ID ,R(40> ,P ( 40) , AA( 40 I 
DIMENSION SZ(ll) 

COMMON/ VALUES/ SI 11) , HI l i ) • B V 111 > * SSI 2 1) , Z , NUM 
COMMON/QUANIT/ZETAI 11) ,CU0,11) ,ZL 
EQUIVALENCE ( B , BV ) , ( BB , C ) 

NUMBER OF MASSES (N) = NUM 
SCALE FACTOR * Z 
MASS CONSTANT (M) * W 
DAMPING CONSTANT (B) = B 
SPRING CONSTANT (K) * S 

1 FORMAT (12) 

2 FORMAT (2X*10F12*6) 

3 FORMAT (2X,2F20>8) 

4 RE AD( 5 , 1 ) ( NUM) 

5 ZL*NUM 

6 LA*2*NUM 

7 M0»2*LA 

8 MA*NUM— I 

9 DO 11 1*1,11 
10 S( U = L.O 

n wm*i.o 

12 DO 13 LUD* 1 , NUM 

13 READ(5,3)(W( LUD ) , S( LUD ) ) 

14 WRITE! 6,3) ( W, S ) 

15 DO 17 1*1, MA 

16 IJK=I+NUM 

17 001 1 JK ) *S( I ) **2 

18 ZW*W(l)*W(2)*W(3)*W(4)*W(5)*W(6)*Wm*W(8)*W(9)*W( 10) 

19 ZS *S I 1 )*Sl2>*SI3)*S(4)*S(5)+S(6)/ZW*S<7i*S(8)*S<91*S(10> 

20 DO 23 1*1,100 

21 ZS*ZS/10. 

22 IFCZS-10* >24,24,23 

23 CONTINUE 

24 PL= I 

25 Z*— PL/ ZL 

26 WRITE! 6, 3) ( ZW,PL»ZL»Z, ZS) 

27 DO 29 1*1,10 

28 V( I)=1.0 

29 ZETAI I ) *1.0 

30 DO 32 IA*1 , 10 

31 DO 32 JA=1 , 10 

32 BB ( I A, JA ) *0*0 

33 DO 63 I V Y* 1 , NUM 

34 DO 62 IZZ*1 , NUM 

35 DO 36 IM*1,10 

36 B( IW)*0.0 

37 B( IZZ)*ZETA( IZZ) 

38 CALL COEFF 

39 DO 40 LP=1 » 1 1 

40 SZ(LP) =SS( 2*LP-1> 

41 CALL POLYRT tSZ,R, NUM , 1 ) 

42 DO 43 I V*1 , LA 

43 R( I V ) = SORT ( ABS t R( I V ) )) * 10 . ** ( -Z/ 2 . ) 

44 MX=2*MA 

45 DO 51 K=2 , MX , 2 

46 DU 51 J=K,MX,2 

47 IF(ABS(R(K) )-A8S( R( J+2 ) ) )51, 51,48 

48 P(K)=R(J+2) 

49 R( J + 2 ) =R t K ) 

50 R( K ) *P ( K ) 

51 CONTINUE 

52 DO 53 LL* 1 , MA 

53 00<LL)*IB(LU*R(2*IVY) )**2 

54 UY = 1 • 

55 DO 57 1*1, MA 

56 I DZ= I+NUM 

57 UY*SURT( 00( I ) +0U( IDZ ) ) *0Y* 10 • ** ( ( 2 . *ZL- 1 . ) *Z / ( 2 . * { ZL-1. ) ) ) 

58 0M=R(2*IVY)*10.**( Z/2. ) 

59 W0*0M**2 

60 BEQ=( l./OY) *( ( ( ( SS( 10) *WO-SS( 8) )*WO+SS( 6) J*W0-SS(4) )*W0+SS12I ) 

61 BEQ=ZW*ABS(BEQ) 

62 BB( IVY, IZZ)*BtQ/B( III) 

63 CONTINUE 

64 WRITE(6,2) (BB) 

65 CALL GAUSSR 

66 STOP 

67 ENO 
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DEVELOPMENT OF METHOD TO OBTAIN COEFFICIENTS OF 
CHARACTERISTIC POLYNOMIAL 


The method used is a modification of Leverriers method, as found in reference 3. 
Equation (5) may be rewritten in matrix -vector form such that 


[[X > 2 + [B]S + [K] [x(S)] = [F(S)] 


(Dl) 


Then, to obtain the characteristic equation, [F(S)] is set equal to zero, and equation 
(Dl) is reduced with the substitutions 


to the form 


<Pj(S) = x(S) 
x( S) = S^ 1 (S) = cp 2 ( S) 


|>(S)] = [A][^(S)] 


(D2) 


(D3) 


(D4) 


where 


[*«] = 


<Pl(S) 

<p 2 (S) 


(D5) 


[A] = 


[O] 


[I] 


-MM]- 1 -[HIM]- 1 


(D6) 


The coefficients are then found as a function of the trace of the powers of the [A] 
matrix. Symbolically, the general term is 


2n-k 


p k ~ ~ £ (Pk-l?l + p k-2 r 2 + • • • + Pl?k-1 + r k) k - 1, 2, . . . , 2n (D7) 

y £ = trace [A]^ (D8) 
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Then 


a 2n-l = ~y\ 

(09) 

a 2n-2 = ‘ ^ (p l^l + y?) 

(DIO) 

2n-3 = " ^ ( p 2 r l + p l y 2 + y $ 

(Dll) 

a 2n-k = Pk 

(012) 

a 2n = 1 

(013) 


The 2n+l coefficients of the characteristic polynomial may be generated quite rapidly by 
this repetitive procedure. 

The flow charts for COEFF and GAUSSR are given in figures 10 and 11, respectively. 
The subroutines are presented after figures 10 and 11. 



Figure 10. - Flow chart for subroutine COEFF. 



Figure 11. - Flow chart for subroutine GAUSSR. 
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SUBROUTINE COtFF 

DOUBLE PRECISION A , B, C , S IGMA , CC, CS, P * t3 Z , SZ , WZ 

DIMENSION AC2u,2C ) ,8120,20) , CC 20 , 20 ) , CS l 2 1 ) , CC C 20 ) , SIGMA (20) , 
lU(ll),BZ(ll>,SZ(ll),WZ(ll> 

COMMON /VALUtS/$( l 1 ) , W ( 1 1 ) , BV t 1 1 ) , SS ( 2 1 ) , Z,NJM 
C DAMPING CONSTANT (B) - 6V 

C SCALED MASS CONSTANT = WZ 

C SCALEO DAMPING CONSTANT = 52 

C SCALED SPRING CONSTANT = SZ 

C RHO(L)=CC(L) 

1 LU-NUM- l 

2 LY *NUM 

3 LW=2*LY 

4 LX=NUM+l 

5 DO 8 1= i » L Y 

6 BZ(I + l)«(10.**<Z/2.n«BV(It 

7 SZ(l+l) = (iO.**Z)*Sm 

8 WZl 1*1 )=W( I ) 

9 BZ(i)=Q.OOQ 

10 SZ(1)=O.ODO 

11 $suw+n = i.o 

12 CS(1)=1.0 

13 00 15 L*1,LW 

14 CC ( L ) =0 • COO 

15 SIGMAC L)=Q.ODO 

16 DO 19 1=1, LW 

17 DO 19 4= 1 , L W 

18 AC 1,41 =0.000 

19 C( I , 4 ) =0*000 

20 CONTINUE 

21 DO 26 1=1, LY 

22 LV=I+LY 

23 AC I, LV 1 = 1*0 

24 AC LV » L V ) =- C BZ ( 1 1 +BZ( 1+1 ) ) / WZ ( I+l ) 

2 5 ACLV,I)=-CSZm+SZCl + lM/WZ( I + l) 

26 CONTINUE 

27 DO 33 1*1, LU 

28 LV=I+LY 

29 A(LV*l,LV)*BZ( I + D/WZU + i) 

30 ACLV,LV+1)=BZ( I+1)/WZ( U2I 

31 ACLV,i+l)=SZ( I+1)/WZ( 1+2) 

32 ACLV + 1, I > ~ S Z C I+l ) /WZ( 1 + 1) 

33 CONTINUE 

34 DO 36 1=1, LW 

35 DO 36 4=1, LW 

36 B( 1,4) =A( I, J) 

37 DO 38 1*1, LW 

38 S1GMAC 1)=AC I, D+SIGMAC l) 

39 CC C 1 ) *— S IGMA ( l ) 

40 SS ( LW) *CCC 1 ) 

41 DO 71 L=2 , LW 

42 DO 59 K= 1 * LW 

43 DO 45 1*1, LY 

44 LV=I+LY 

45 CC I,K)=AC I,LV)*BCLV,K)+C< 1,K) 

46 DO 49 1=LX , LW 

47 LV= I-L Y 

48 CC I,K)*AC I,LV)*BCLV,K)+C( I,K) 

49 CC I , K) =A ( 1 , I)*BC1,K)+C( I,K) 

50 LR =L W- 1 

51 DO 54 1=LX,LR 

52 LP=I-LU 

53 CC I,K)=A(1*LP)*B(LP»K)+CCI,K) 

54 CC 1,K)*AC l,l+i)*BC I+1,KI+C< 1,K) 

55 L0*2+LY 

56 DO 59 I*LO, LW 

57 LC= I-L X 

58 CC I * K ) *A( I,LC)*BCLC,K)+CC I , K ) 

59 CC I,K)*ACI,1-1)*BCI~1,K)+CCI,K) 

60 DO 64 1*1, LW 

61 SIGMACL>=C( 1, D+SIGMACL) 

62 DO 64 4= 1 , L W 

63 Bt 1 , 4 ) =C C 1,4) 

64 CC I , 4 ) =0*0 DO 

65 P*L 

66 DO 71 J= 1 ,L 

67 N=L~J+1 

68 M*LW+1~L 

69 CS C L ) *CC C L“ 1 ) 

70 CCC L )*-( 1*D0/P)»CCS(N) *S1GMAC4 ) ) +CCC L ) 

71 SSCM)=CCCL) 

73 FORMAT (2X,4E20*7) 

74 RETURN 

75 END 


) 

1 

( 
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SUBROUTINE GAUSsR 

DIMENSION AUO,ll>,aill> .BB(IG) 

CUMMUN/QUAN I T/ZE TA( 11 ) ,Ct 10* 11) * ZL 
LQUI VALENCE t B , 2ETA) 

C C M A TR I X= MATRIX OF INFLUENCE COEFFICIENTS 

C READ STATEMENT AUGUMENTS C MATRIX, READS B EQUIVALENT 

1 FORMAT ( 2X » 5F 1 5 • 8 ) 

2 FORMAT (2X,5E22.B) 

3 NUM=ZL 

4 LX =NUM+ 1 

5 READ(5,l)(Cll»LX)fCt2*LX)fC(3fLX)tCI4,LX),C(5»LXJfCt6»LX)fC(7» L*X) 
61,C(8,LX) ,C(9,LX) ,C<10,LX>) 

7 WRITEt6,2) tC,NUM) 

8 LU=NUM-L 

9 DO 11 1=1, NUM 

10 DO 11 J= 1 , LX 

11 A ( I, J)=C( I,J) 

12 DO 38 L=i,LU 

13 DO 14 I I =L , NUM 

14 Qd(lII=MII,L) 

15 DO 21 KK=L , LU 

16 DO 21 J J=KK , LU 

17 IFCA6S(BB(KKII-ABS{BBC JJ+1) ) ) 21,21* 18 

18 2= BB t KK ) 

19 BB(KK>=BB< JJ+1) 

20 BBIJJ+1)=2 

21 CONTINUE 

22 00 25 JK*1 * NUM 

23 IF(ABS(A(JK»L) )— ABS( BB ( NUM) > >25,24,25 

24 L2 = JK 

25 CONTINUE 

26 DO 30 I H=1 * LX 

27 2= AI L2, IH) 

28 A { L Z , I H ) =A ( L , IH) 

29 A(L, I H ) =2 

30 CONTINUE 

31 DO 38 1=1, LX 

32 B ( I)=0.0 

33 K= NUM+ 2- I 

34 DO 35 J=L,NUM 

35 At J,K)*A(J,K)/A(J,L) 

36 DO 37 J=L , LU 

37 At J+1,K)=A( J+1,K)-A( L,K) 

38 CONTINUE 

39 DO 41 M= 1 , NUM 

40 I =LX-M 

41 Bt I ) =A ( 1 ,LX )/A(I, 1 1 -At I f 2)*B(2)-At I, 3)*B(3)-AC I ,4)*Bt4)-A( 1,51* 
421 B( 5) -A l 1 ,6) *B( 6) -At I ,7) *B(7)-A( I ,8) *B( 8) -At I ,9) *Bt 9) -At I, 10) *Bt 10) 

43 WRlTEt 6,44) t B) 

44 FORMAT ( 2X,E20*8) 

RETURN 

END 
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